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The theory to reconstruct the atomistic-level chain diffusion from the accelerated dynamics that is measured 
in mesoscale simulations of the coarse-grained system, is applied here to the dynamics of ds-l,4-Polybutadiene 
melts where each chain is described as a soft interacting colloidal particle. The rescaling formalism accounts for 
the corrections in the dynamics due to the change in entropy and the change in friction that are a consequence 
of the coarse-graining procedure. By including these two corrections the dynamics is rescaled to reproduce the 
realistic dynamics of the system described at the atomistic level. The rescaled diffusion coefficient obtained 
from mesoscale simulations of coarse-grained cis-l,4-Polybutadiene melts shows good agreement with data 
from united atom simulations performed by Tsolou et al. The derived monomer friction coefficient is used as 
an input to the theory for cooperative dynamics that describes the internal dynamics of a polymer moving 
in a transient regions of slow cooperative motion in a liquid of macromolecules. Theoretically predicted time 
correlation functions show good agreement with simulations in the whole range of length and time scales in 
which data are available. The theory provides, from data of mesoscale simulations of soft spheres, the correct 
atomistic- level dynamics, having as solo input static quantities. 



I. INTRODUCTION 

Computer simulations of macromolecular liquids de- 
scribed at the atomistic level are extremely useful 
because they bridge information between the micro- 
scopic molecular structure of the polymeric system, at 
the given thermodynamic conditions, and its macro- 
scopic properties, such as viscosity, diffusion, and 
dynamical- mechanical response, which are observed 
experimentally^ 2 It is unfortunate that, despite the 
progress already occurred in the computational hardware 
and software, the study of the dynamics of polymeric liq- 
uids by atomistic computer simulations is still limited by 
the impossibility of simulating polymer dynamics in the 
wide range of time and length scales relevant for these 
systems. It is known that long simulation trajectories 
deteriorate with time, following a law defined by the Lya- 
punov exponent characteristic of the system, and the re- 
sulting long-time dynamics is affected by errors^ 

Atomistic simulations of polymeric liquids are limited 
either in the length and number of macromolecules, or in 
the maximum timescale that can be reached. The long 
timescale is a regime of great interest for these systems, 
because of the relevance of dynamical-mechanical prop- 
erties of macromolecules in that regime for their indus- 
trial applications. Because the limitation concerns only 
the region of long-time, large-scale properties, it is pos- 
sible to shift the focus of a simulation to this region of 
interest by the reduction of the simulated degrees of free- 
dom through the averaging of local scale properties, i.e. 
adopting a so-called "coarse-graining" procedure. 

While simulations of coarse-grained systems afford 
to represent larger length- and longer time-scales than 
atomistic simulations, they also are marred by the short- 
coming that the measured dynamics is unphysically fast. 
This is the advantage that allows one to perform simu- 



lations in the long-time regime, but it is also a problem 
because the dynamics in the simulation becomes too fast 
and measured diffusion coefficients are too large. De- 
pending on the level of coarse-graining the dynamics can 
become orders of magnitude faster than the atomistic 
dynamics 

Being aware that the dynamics has to be rescaled 
to recover the correct timescale^— one resorts to the 
most conventional method, which is to build a calibration 
curve calculated by superimposing the long-time dynam- 
ics of the coarse-grained and atomistic simulations ! 1Q T n 
The hope is for this calibration curve to be transferable 
to other systems, or to be identical, and so applicable, 
for thermodynamic conditions close to, but different than 
the phase point for which the curve was optimized in the 
first place. Otherwise, if an atomistic simulation has to 
be ran every time we need to rescale a coarse-grained 
simulation there is no advantage in running simulations 
of coarse-grained (CG) systems. The rescaling factor of 
the conventional calibration curve is, however, purely a 
numerical correction. There is no physical motivation for 
one to assume that this correction has to be identical for 
any system or thermodynamic condition different from 
the ones in which it was optimized. Given that the ef- 
fective potentials between coarse-grained units are free 
energies, and as such they are parameter dependent, it is 
very unlikely that numerically optimized corrections will 
be identical for different systems. However, if the calibra- 
tion curve is applied for conditions of temperature, den- 
sity, and degree of polymerization, very close in the phase 
diagram to the ones in which the calibration curve was 
originally optimized, there is hope that the error is very 
contained. Alternatively, in the case that very limited 
levels of coarse-graining are applied, e.g. combination 
of a few atoms to obtain a new coarse-grained unit, the 
error incurred in using the calibration curve in different 
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conditions, can be small enough to be within the numer- 
ical error of the calculation. This is in fact the reason for 
the observed good agreement in the rescaled dynamics of 
systems with contained level of coarse-graining*^ — 

The goal of this paper is orthogonal to the more con- 
ventional calculation of an optimized calibration curve 
for systems with limited level of coarse-graining. In fact 
our goal is to investigate the physical motivation that 
leads to the accelerated dynamics in first place and to 
provide, through our analysis, a theoretical approach to 
calculate, from first principle equations, the correction 
factor that needs to be applied to the fast dynamics mea- 
sured in coarse-grained simulations to recover the slower 
and more realistic dynamics of the atomistic description. 
We study the most extreme level of molecular coarse- 
graining, where one polymer chain is represented as a 
soft sphere, because i) the description is analytical al- 
lowing the formal study of the rescaling problem, ii) it 
provides a solid test of the reconstruction procedure be- 
cause possible errors would be maximized, iii) this level 
of coarse-graining is associated with the largest compu- 
tational gain, so it is important. 

Recently we have proposed an original procedure to 
rescale the dynamics measured in simulations of coarse- 
grained polymer melts AS The procedure we proposed 
does not require to run atomistic simulations to calibrate 
the dynamics of the coarse-grained simulation. Once 
one parameter, specifically the diameter of the hard-core 
monomer potential, is fixed, the rescaling factors are fully 
determined. The advantage of this method is clear, given 
that the coarse-grained simulation can reach larger-scale 
longer-time properties than the atomistic simulation. Be- 
cause no atomistic simulations are needed, the measured 
dynamics can be directly rescaled to obtain the dynamics 
of the real system we wanted to study. 

The method has been formally derived, and then ap- 
plied to systems of polyethylene (PE) melts. The choice 
of PE as a test system was motivated by the wealth of 
experimental and simulation data available in the lit- 
erature. We have shown that the proposed theory is 
able to predict diffusion coefficients in agreement with 
both atomistic simulations and experiments for systems 
in both the unentangled and entangled regimes, for a 
range of temperatures and densities. For entangled sys- 
tem the approach of dynamical reconstruction includes 
an extra loop in the calculation of the friction coefficient, 
which accounts for the fact that in real systems described 
at the atomistic level, the dynamics is slowed down by 
the presence of entanglements. We hypothesized for the 
entanglements to relax in time following the same dy- 
namical mechanism of single chain inter diffusion. In this 
paper entangled chain dynamics is not investigated as 
most of the samples in the atomistic simulations are ei- 
ther unentangled or in the region of crossover from the 
unentangled to the entangled regimes. 

One of the questions that still needed to be answered 
was if our method was able to generate high quality 
predictions also for macromolecular liquids with a more 



complex monomeric structure than the simple polyethy- 
lene. A first study of this issue is presented in this 
manuscript which investigates the agreement between 
theoretically predicted diffusion coefficients for ds-1,4- 
Polybutadiene (PB) chains and data from molecular dy- 
namics simulations performed by Tsolou, Mavrantzas 
and coworkers J^— The approach is completely general 
and can be applied to polybutadiene chains with differ- 
ent tacticity. For these systems the difference in local 
chain semiflexibility, which corresponds to different over- 
all chain dimension, i.e. a different radius-of-gyration 
(R g ) and statistical segment length (/), ultimately leads 
to a different diffusion coefficients for polybutadienes 
with different tacticity. In this paper we limit our cal- 
culations to cis-l,4-Polybutadiene samples because those 
are the samples investigated by atomistic simulations and 
provide a effective test of our procedure for dynamical re- 
construction. 

Our rescaling method considers two contributions 
emerging from the consequence of applying the coarse- 
graining process and derives the theoretical corrections 
that need to be applied to the dynamics of the coarse- 
grained description to recover the modalities of the orig- 
inal atomistic system. To calculate the analytical cor- 
rection we adopt for the atomistic descriptions a simple 
bead-and-spring description of the chain where each PB 
polymer of N monomeric units is represented as a col- 
lection of N beads connected by semi flexible springs of 
length I = yj 6 /NR g , and the overall chain dimension 
is defined by the radius-of-gyration of the molecule, R g . 
This model is a Rouse approach modified to include chain 
semiflexibility and has been shown to represent well the 
diffusion of polymer melts in the unentangled regimeJ^ 
From the comparison of this model and the soft sphere 
representation of the chain, it is possible to evaluate the 
two contributions that are responsible for the accelerated 
dynamics, and which have to be included to correct the 
measured dynamics and bring it to the realistic values of 
the atomistic (here bead-and-spring) description. 

The first correction emerges from the consideration 
that due to the process of averaging the intramolecu- 
lar degrees of freedom, the coarse-grained system experi- 
ences a change in its entropy, because a number of local 
states are neglected. While the atomistic system devotes 
time to explore local configurations, the coarse-grained 
system doesn't. To recover the correct dynamics those 
states need to be included a posteriori in the form of an 
entropic contribution to the free energy that rescales the 
time measured in the CG simulations. 

The second correction comes from the change of the 
"shape" of the molecule once it is coarse-grained. Macro- 
molecules described at the atomistic level become rep- 
resented by chains of effective atoms, or as bead-and- 
springs, chains of soft blobs, or even each molecule as a 
soft sphere. In all these multiple forms the surface of each 
molecule available to the surrounding molecules, i.e. the 
"solvent", is different. The hydrodynamic radius, rhydr-, 
of the coarse-grained unit in each description is differ- 
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ent and so is its friction coefficient, defined by Stokes' 
law as £ = 67rr] s rh y dr, where r] s is the viscosity of the 
medium. By solving the memory function for the friction 
coefficient in each description we are able to calculate the 
correction factor that has to be applied to each friction 
coefficient to recover the description at a different level 
of coarse-graining. 

In this paper our rescaling method is first briefly sum- 
marized and then applied to study the dynamics of cis- 
1,4-PB and compared with simulation data by Tsolou 
et.al.— who also investigated the effect of temperature 
and pressure on the simulations ] 16 i 17 The global dynam- 
ics of the polymer is represented by the diffusion coeffi- 
cient and the decay of the self-correlation for the end-to- 
end molecular vector. These dynamical quantities are re- 
constructed using the method by Lyubimov at al.-r^ from 
the information contained in the trajectories of the meso- 
cale simulations of the coarse-grained polymer melts and 
show good agreement with the data from united atom 
simulations. For the internal dynamics the theory for 
cooperative motion^"— is used to predict monomer dy- 
namics and the dynamical structure factors, which com- 
pare well with the simulations. If atomistic simulations 
are not available, the theory of dynamical reconstruction 
is able to predict the correct dynamics from the rescal- 
ing of the mesoscale simulations having as an input the 
polymer radius-of-gyration, which can be calculated from 
simple structural models for polymers. The accuracy of 
adopting a freely rotating chain model for the calculation 
of the polymer radius of gyration is discussed in the last 
section. A brief summary concludes the paper. 



II. RESCALING OF THE FREE-ENERGY AND 
RESCALING OF THE FRICTION 

In this section we briefly present an overview of the 
main steps in our procedure. A complete and detailed 
presentation of the same has been already published^ 
and will not be repeated here. We start from the con- 
sideration that the procedure of coarse-graining corre- 
sponds to eliminating some internal degrees of freedom, 
by combining groups of atoms into one effective CG unit. 
Specifically, our rescaling theory has been developed for 
an extended level of coarse-graining where intramolecu- 
lar coordinates are fully averaged out and the polymer is 
represented as an isotropic sphere centered on the center- 
of-mass (com) of the macromolecule. There are some ad- 
vantages in the choice of this description. First, the CG 
representation is fully isotropic, and even if it is known 
that the shape of a polymer is closer to an ellipsoid than 
to a sphere^ the total correlation function of the poly- 
meric liquid, i.e. the structure of the polymeric liquid, 
statistically averaged over all the possible configurations 
is well reproduced by this model; second the formalism 
is analytical, which allows us to derive formally many of 
the physical quantities of interest for both the static and 
the dynamic properties of the coarse-grained system as 



a function of the atomistic description; third all the con- 
sequences of the coarse-graining procedure are enhanced 
with this extreme level of coarse-graining and are easier 
to study in such a description. If the procedure adopted 
were not precise, the error would be maximized because 
of the high level of coarse-gaining; the fact that the the- 
ory, which is predictive, is found to be in quantitative 
agreement with both simulated and experimental data is 
encouraging^ 

One disadvantage of having such a coarse-grained rep- 
resentation of the molecule is that no information is col- 
lected from the mesoscale simulation on the internal dy- 
namics of the polymer. In this way, the soft sphere 
representation only provides information on the diffu- 
sion coefficient of the molecule. However, we show in 
this paper that from the knowledge of the diffusion co- 
efficient and the related monomer friction coefficient for 
unentangled chains it is possible to recover correctly the 
internal dynamics of the chain by applying our theory 
for the cooperative dynamics of a group of interacting 
chains, which calculates the single chain dynamics for a 
macromolecule whose dynamics is coupled by the pres- 
ence of intermolecular forces with the other interpene- 
trating macromolecules in the liquid. 

In every CG model, the averaging of the intramolecular 
degrees of freedom leads to a speeding up of the dynam- 
ics. When a polymer is represented as a soft sphere the 
dynamics in the mesoscopic simulation is orders of mag- 
nitude faster than in the atomistic description. One of 
the reasons for this acceleration is the fact that while a 
macromolecule needs to explore a large number of inter- 
nal chain configurations for each position of the center-of- 
mass, the soft sphere instead is free to move in the three- 
dimensional space undergoing Brownian motion j 23 i 24 To 
take into account the missing contribution of the time 
necessary to explore the internal degrees of freedom, the 
time measured in the mesoscale simulation of the soft 
sphere is multiplied by a rescaling factor that is calcu- 
lated as the ratio of the energy due to the internal de- 
grees of freedom in the two representations, i.e an atom- 
istic and a soft sphere representations. For the atomistic 
representation we adopt an analytical bead-and-spring 
model, as previously described^ This model has been 
shown to provide a quantitative description of the dy- 
namics for polymeric liquids and for proteins in theta 
solvent £L£L 

The first correction term is calculated from the def- 
inition of the internal energy for a simple bead-and- 
spring model and for the soft sphere CG model. The 
mesoscale (MS) molecular dynamic (MD) simulations are 
performed using a potential expressed in dimensionless 
units, which is a common procedure. The time from the 
simulations has to be properly normalized so that the di- 
mensionless simulation time from MS simulations, once 
dimensionalized and rescaled reads as 



where m is the molecular mass of one chain and R g is 
the radius-of-gyration of the molecule, which is an in- 
put quantity in our approach. The internal degree-of- 
freedom averaged out in the coarse-grained description 
are accounted for by introducing the |iV factor where N 
is the number of beads. 

The second rescaling is due to the change in the shape 
of the molecule represented in the two levels of coarse- 
graining, and so the effective friction in the two repre- 
sentations. In the long time the com mean-square dis- 
placement is related to the molecular diffusion coefficient 
as 



The contributions that are still unknown are the friction 
coefficients in the atomistic and coarse-grained descrip- 
tions, which can be solved analytically starting from the 
definition of the memory function^ 

For a bead-and-spring model the monomer (bead) fric- 
tion is given by: 

-j N poo 

^ - TV E / d ^a 9 b(r) . (4) 

a, 6=1 J ° 

where the memory function is defined as2£ 



{AR 2 (t)) =6Dt , (2) 

where D = ksT '/ '(N£ m ) is the diffusion coefficient of the 
system, which has to be calculated. If we use as a starting 
point the diffusion coefficient measured in the mesoscale 
simulation D MS = k B T/( soft , the mean square displace- 
ment in the rescaled formalism is give by 



(AR 2 (t)) =6D MS j^t , (3) 

where the problem reduces to the calculation of the fric- 
tion coefficient in the monomer and soft-sphere represen- 
tations. The friction coefficients are calculated from the 
solution of the memory function in the two representa- 
tions (bead-and-spring and soft-sphere). 

The formalism presented so far, adopts a Rouse-like 
description of the single chain diffusion for the atomistic 
level representation of the system. At the coarse-grained 
level the chain is represented as a soft-sphere. Eqj3] shows 
how the diffusion coefficient measured in a mesoscale sim- 
ulation of a liquid of soft-spheres, D MS , needs to be 
rescaled to give the mean-square displacement of the real 
chain, (AR 2 (t)}, as a function of the rescaled time, t. 



r o ,6(t) = yj drj dr'g(r)g(r')F(r)F(r')r-r' 

x [dRS% b (R;t)S Q (\r-r' + H\;t), (5) 

where f3 = P is the monomer density and f, ?' 

are the unit vectors characterizing the direction of the 
forces acting on monomers a and 6, respectively. S®(r, t) 
is the projected dynamic structure factor, which includes 
both intra- and inter-molecular contributions, i.e. inco- 
herent and coherent scattering, and is approximated by 
its non-projected form as S®(r,t) ~ S(r,i), which is an 
acceptable approximation when the Langevin equation is 
expressed as a function of the slow variables. 6 

In Eq.(|5]) g(r) is the monomer pair distribution func- 
tion of the molecular liquid, and F(r) is the force due 
to the effective potential between two monomers, ob- 
tained from the solution of the Ornstein-Zernike (OZ) 
equation by applying the Percus Yevick closure. The 
monomer potential, which in the atomistic-level simula- 
tion is a Lennard- Jones potential, is approximated here 
by an effective hard-sphere with an effective diameter, 
d, to mimic the properties of the L-J potential. All the 
physical quantities that appear in the equation are known 
and the analytical expression for the monomer friction is£ 
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This equation contains quantities that are well defined 
once the system of interest is selected. For example 
g(d) is the pair distribution function at contact, £ p = 
R g /(y/2 + 27T Rgp/N) is the density fluctuation correla- 
tion length with p, and N already defined in the text, 
and ho is defined as ho = h(k = 0) = — (1 — 2^/ Rg) / p c h- 
The only physical quantity that needs to be determined 
is the hard-sphere diameter, d. This is defined once the 
Lennard- Jones potential is mapped onto an effective re- 
pulsive hard-sphere system as described later in this sec- 
tion of the paper. 

In the soft colloidal representation, the friction coeffi- 
cient is analytically calculated from the definition of the 
memory function as 

Csoft = (P/3)Pch j~ dt f dv J dv f g(r)g(r f )F(r)F(T f )v-v' 

x J dHS(R-t)S(\r-r f + R|;t), (7) 

where p c h = p/N is the chain density. All the other 
physical quantities that appear in the integral, for exam- 
ple the pair distribution function, g(r), and its related 
force, F(r), are now defined in relation to the descrip- 
tion of the polymeric liquid as a liquid of soft spheres, 
which are point particles interacting through a soft re- 
pulsive potential. Each of these quantities is analytical 
and calculated from the solution of the OZ equation with 
the hypernetted chain closure approximation^ 

The regime of interest in our calculations is the dif- 
fusive limit where, in reciprocal space, the wave vector 
k <C 1/Rg- Here the dynamic structure factor can be 
approximated as 
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with the total correlation function of the soft-sphere rep- 
resentation is given in the limit of long chains (N > 30) 
by the approximated expression 27 



h soft(r) 



The resulting equation for soft particle friction is£ 
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where the value of each physical parameter (temperature 
and density) and molecular parameters (degree of poly- 
merization and radius-of-gyration) that enter this equa- 
tion is defined once the system to simulate is selected. 
These parameters are identical to the ones that are input 
to the equation of the monomer friction, Eq.©, given 
that the two equations are representations for different 
levels of coarse-graining of the same system. 

When the two equations for the monomer and the soft 
colloidal friction coefficients are introduced in Eq.(|3]) the 
diffusion coefficient for a single chain in a liquid is recov- 
ered as 



D = D MS X soft /(N( m R g ^3mN/(2k B T)) , (11) 

where we have assumed that the diffusion coefficients 
that appear in Eq|6]and EqJTD]are identical, as the long- 
time relaxation of the dynamic structure factor in both 
CG descriptions is guided by molecular diffusion. 

To summarize, as the first step the thermodynamic 
and molecular structure, i.e. the radius of gyration or 
equivalently the persistence length, of a polymer melt has 
to be defined. Then, MS MD simulations are performed 
for a liquid of point particles interacting through a soft 
pair potential, (3v so f t (r), as described in several of our 
previous papers, where 

Pv 8oft (r) = h 8oft (r) - ln[l + h soft (r)} - c soft (r) , (12) 

with 



Csoft(k) = h ao f t (k)[l + pchhsoft(k)]' 



(13) 



Rq 



The mean-square-displacement of the soft spheres is mea- 
sured from the MS MD trajectories as a function of time 
(dimensionless), and the resulting diffusion coefficient is 
rescaled following EqJTT]to obtain the reconstructed dif- 
fusion coefficient. 

For the solution of the rescaling factors in EqJTT] the 
only parameter that has to be defined is the value of 
the effective hard sphere diameter, d, in Eqj6j In our 
formalism the elementary interaction between monomers 
belonging to different chains, which in the UA MD is 
a Lennard- Jones potential, is approximated by a hard- 
sphere interaction with an effective diameter, d, to allow 
for the analytical solution of the monomer memory func- 
tion. Under fixed thermodynamic conditions, the value 
of d should depend only on the local monomer structure, 
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d [A] d [A] 

FIG. 1. Dimensionless monomer friction coefficient as a func- 
tion of the hard sphere diameter, based on Eq.JHj). a) cis- 
1,4-Polybutadiene samples with N= 32, 56, 128, 320 (solid, 
dashed, dot-dashed, and dotted lines correspondingly) and 
parameters as reported in Tableland ref.— . b) Polyethylene 
samples with N=30, 44, 96, 270 (solid, dashed, dot-dashed, 
and dotted lines correspondingly) and parameters as reported 
in refs.— and^. Horizontal lines represent 1/N values, following 
the diffusion coefficient for unentangled chains, Df3( m = 1/N. 



and be independent of the degree of polymerization, be- 
cause the monomer interaction potential is a local prop- 
erty of the chain. Therefore in our model, once the value 
of d is chosen, this value is kept fixed for all samples with 
different molecular weights. Note that for convenience 
we use the term monomer to identify the CH X group, 
where x = 1, 2 or 3. 

To evaluate the value of d we analyze the dimension- 
less quantity D( m /(kBT), which is identical to TV -1 if 
the chain obeys strictly Rouse dynamics in the long- 
time limit. To fulfill this requirement the chain has to 
have N much smaller than the entanglement degree-of- 
polymerization N e and larger than the value of N = 30 
for which chains start to obey Gaussian statistics, i.e. the 
central limit theorem applies. 

Figure [T] a) displays the dimensionless friction coeffi- 
cient, Df3( m , from Eq.((6j) as a function of the hard sphere 
diameter, <i, for PB samples of increasing degree of poly- 
merization, N. The parameter values (p m , T, TV, R g ) are 
chosen to match those of UA MD simulations in ref as 
reported in Table HI Horizontal lines represent the 1/N 
value for each sample, which is the Rouse value of the 
physical quantity Df3( m = 1/N in the diffusive regime, 
and reproduces the scaling behavior of the diffusion for 
unentangled, short-chains. For comparison Figure [T] b) 
displays the analogous plot for the PE samples, which is 
reproduced from refs£ and£. 

The change of d with TV observed in Figure [1] is a con- 
sequence of two different effects: for short chains, end ef- 
fects are important, while for long chains the crossover to 
entangled dynamics starts to be important. In this way, 



it could be argued that in our method the actual fitting 
parameter is not but the length of the chosen sam- 
ple for which d is determined. Clearly when entangled 
chains are forced to behave according to the Rouse ex- 
pression the hard sphere diameter needs to be artificially 
decreased in order to compensate the overlooked increase 
in the monomer friction, which is associated with entan- 
glement effects growing with increasing N. Whereas the 
range of lengths in the unentangled region can be rather 
wide, it is always possible to choose the chain length 
knowing experimentally measured or theoretically esti- 
mated value of Ag, 28 and by optimizing our choice of 
d by testing the predictions of the theory against data 
from simulations of short chains. However comparison 
with simulations is not necessary, as the only needed in- 
formation is the value of N e and the value of R g for a 
sample of short, unentangled, chains (N < N e ). 6 

In our previous study of PE melts we fixed the hard 
sphere diameter following the same procedure described 
here, and we selected a sample with N = 44 to cal- 
culate d. The obtained value of d = 2.1 A (see Fig- 
ure [TJd), was larger than the carbon-carbon bond length 
lb = 1.54A, but smaller than the Lennard- Jones param- 
eter a ~ 3. 9 A typically used in UA MD simulations, 
which seems to be a reasonable choice of the parame- 
ter. Compared to PE, PB chains are more flexible and 
so they are less entangled. The estimated entanglement 
degree-of-polymerization, N e , is calculated from the gen- 
eralized formula for the plateau value in the shear relax- 
ation modulus, 28 which gives N e ( shear ^ ~ 330 backbone 
carbon atoms for PB, compared to N e ^ shear ^ ~ 55 for PE. 
Notice that the value of N e from shear measurements is 
different from the value obtained from the analysis of 
scattering experiments which gives, for example, for PE 
samples N e ~ 130. 

It is important to notice that even if there is some 
freedom in choosing the reference sample, the deviations 
of d from the chosen value is not large. However, small 
deviations of the value of d can lead to different values 
of the diffusion coefficient. Figure [2] displays a compari- 
son of the predicted diffusion coefficient as a function of 
N when different values are chosen for the hard-sphere 
parameter d. It is possible to see that for values of d 
obtained by enforcing Rouse dynamics for three different 
chain lengths that obey diffusive Rouse dynamics in the 
long-time regime, the predicted diffusion coefficients are 
all in reasonable agreement with the values of the UA 
MD simulation. 

Notwithstanding the fact that for PB there is a wide 
range of degrees-of-polymerization, 30 < N < 330, for 
which the hard-sphere diameter, could be optimized, 
the best agreement between theory and simulations is ob- 
tained for values of the parameter d optimized for short 
chains, N w 50. This can be explained by considering 
that in PB samples the crossover region from unentan- 
gled to entangled dynamics is extended (see for example 
Figure[3l where deviations from the A --1 scaling occur al- 
ready for N « 100). Samples with A" approaching N e are 
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FIG. 2. Diffusion coefficient predicted from the rescaled MS 
MD, when different values of the hard-sphere diameter are 
chosen as an input to the reconstruction procedure. Assum- 
ing different values of d, calculated by enforcing Rouse dif- 
fusive behavior for N=56 (d=1.4672), 64 (d=1.4051), or 80 
(d= 1.374) (circles, diamond, triangles correspondingly) leads 
to diffusion coefficients in good agreement with the UA MD 
simulation data (filled squares). 



not following the unent angled scaling behavior, which is 
necessary condition for the optimization of d using Rouse- 
like diffusion. 

Considering the broader range of unentangeled PB 
chains, we select for the calculation of the hard sphere 
diameter the PB sample with N = 56. The obtained 
value of d ~ 1.47A (which is smaller than the value for 
PE of 2.1A given the flexibility of PB) is consistent with 
the fact that in PB chains half of the carbon bonds are 
shorter double bond and Idb = 1.34A. 



III. THEORETICAL PREDICTIONS OF 
CENTER-OF-MASS DIFFUSION 

The solution of Eq.(|3j) gives the diffusion coefficient 
for the center-of-mass of a polymer described at the 
atomistic level, having as an input the diffusion coef- 
ficient calculated from the mesoscale simulation of the 
CG polymer liquid. The purpose of this procedure is to 
have ultimately mesoscale simulations that, once prop- 
erly rescaled, can provide directly the diffusion coeffi- 
cient, without need of running atomistic simulations. To 
test our procedure we first run the mesoscale simula- 
tions, then we rescale the diffusion coefficient, and fi- 
nally we compare the predicted values of D against atom- 
istic simulations and/or experiments. Thermodynamic 
and molecular parameters entering our equations have to 
be consistent with the parameters of the system against 
which we compare our approach. 

We study the dynamics of cis-l,4-Polybutadiene (PB) 
melts of increasing degree of polymerization and com- 



pare our results with the data from the simulations per- 
formed by Tsolou et alJ^ The input parameters to our 
theory are displayed in Table HI including the degree-of- 
polymerization, TV, and the monomer density p. Table U 
also includes the mean-square-radius-of-gyration calcu- 
lated from the united atom simulations (Rg) • From 
this radius-of-gyration the value of the semiflexibility pa- 
rameter, g, is derived, as the value that corresponds to 
chains with the desired overall dimension, i.e. R g . In 
this way the parameter g is not an independent param- 
eter but is determined by the R g value. The parameter 
g does not enter the equation, but is used in the cal- 
culations reported in the last section of this paper. It is 
important to notice that the values for R g could be taken 
from experimental data, and that, in principle, no atom- 
istic simulations of the system under study are necessary 
for our approach. 

The UA MD simulations were performed in the NPT 
ensemble, at the constant pressure, P, and tempera- 
ture, T, reported in the tabled Note that the values 
of R g for systems PB160, PB200, PB240 were corrected 
from 257A 2 to 275 A 2 , from 304 A 2 to 340 A 2 and from 
401 A 2 to 410 A 2 correspondingly, after consultation with 
the authors of ref The correction was motivated by 
the inconsistency observed with the data reported in 
reference^ when the calculation of R g is performed us- 
ing a standard model, such as the Freely- Rotating-Chain, 
described later in in this paper. 



TABLE I. Simulation parameters for 1,4-cis PB chains of in- 
creasing lengths. 



System p [sites/ A 3 



PB 


32 





0352375 


45±5 





6237 


PB 


48 





0363767 


70±5 





6243 


PB 


56 





0367159 


85±7 





6339 


PB 


64 





0369744 


95±10 





6246 


PB 


80 





0373425 


125±10 





6374 


PB 


96 





0375921 


152±16 





6394 


PB 


112 





0377723 


184±15 





6490 


PB 


128 





0379087 


215±18 





6544 


PB 


140 





0379910 


234±18 





6523 


PB 


160 





0381012 


275±25 





6595 


PB 


200 





0382567 


340±20 





6551 


PB 


240 





0383610 


410±30 





6557 


PB 


320 





0384923 


576±30 





6695 


PB 


400 





0385714 


678±30 





6519 


T= 


413K, P=l atm 









In the mesoscale simulations of a polymer melt, each 
chain is represented as a point particle interacting 
through a soft-core potential derived from the solution of 
the Ornstein Zernike equation applying the Hyper Netted 
Chain closure. MS MD simulations were implemented 
in the microcanonical [NVE) ensemble on a cubic box 



8 



with periodic boundary conditions. We used reduced 
units such that all the units of length were scaled by 
R g (r* = r/Rg) and energies were scaled by /c#T. More 
details of our simulation procedure have been reported 
m previous papers i 

Table [XT] reports the diffusion coefficient directly cal- 
culated from the MS MD in the soft sphere represen- 
tation, -D^ s , and the diffusion coefficient calculated us- 
ing the dynamical reconstruction procedure, D. Finally, 
for comparison, the Table shows the diffusion coefficient 
measured in the united atom simulation, D UA . The diffu- 
sion coefficient measured in the mesoscale simulations is 
several orders of magnitude faster than in the atomistic 
simulations. However, once it is rescaled the diffusion 
coefficient becomes very similar to the one obtained di- 
rectly from the atomistic simulation. It is important to 
notice that, once the parameter d, which is characteris- 
tic of the polymer considered, is determined the diffusion 
coefficient is calculated without any input from the dy- 
namics of the atomistic simulations, so the procedure is 
predictive. 

The samples here are relatively short and these cal- 
culations are performed for unentangled and slightly en- 
tangled systems. For strongly entangled systems, we pro- 
posed in a previous paper a perturbative version of our 
approach that accounts for the fact that both the tagged 
chain and the surrounding chains, that provide the en- 
tanglements, relax following the same dynamics £ In this 
manuscript the PB samples are in the unentangled and 
slightly entangled regimes, which are well represented by 
the theory without the one-loop perturbation. 

TABLE II. Diffusion coefficient reconstructed from MS MD 
simulation compared against UA MD simulations. 



N 


A MS [A 2 /ns] 


D[A 2 /ns] 


D UA [k 2 /ns] 


32 


3875 


184.5 


152.2 


48 


3400 


68.1 


69.4 


56 


3425 


51.9 


46.4 


64 


3275 


37.3 


40.1 


80 


3174 


25.1 


24.9 


96 


3114 


18.0 


23.2 


112 


3235 


15.1 


16.7 


128 


3283 


12.6 


15.5 


140 


3221 


10.5 


12.0 


160 


2963 


9.1 


11.0 


200 


3245 


6.0 


6.6 


240 


3012 


4.5 


5.4 


320 


3386 


3.2 


2.6 


400 


3016 


1.8 


1.6 


d = 


1.4672A(PB56) 





On Figure [3] the diffusion coefficient is presented as a 
function of degree of polymerization N. Filled symbols 
represent UA MD simulations from ref. 15 and open sym- 
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FIG. 3. Center of mass self-diffusion coefficient as a func- 
tion of degree of polymerization, N, for cis-l,4-Polybutadiene 
melts with parameters denned in Table|U Diffusion coefficients 
reconstructed from MS MD by applying our procedure (open 
symbol) are compared against UA MD data (filled symbol) 
from reference^. In analogy with the figure from the source, 
three scaling regimes in terms of power law dependence of 
D oc N b are shown as dashed (b > 1), solid (b « 1), and 
dot-dot-dashed (b ~ 2) lines. 



bols represent MS MD simulations rescaled as in Eq.((3j). 
In analogy with the figure in Ref.— from which the UA 
MD data are taken, we report three scaling regimes in 
terms of power law dependence of D oc TV -6 . For N < 80, 
there is a faster than Rouse regime with D with b > 1, 
which is attributed to the free- volume effects due to chain 
ends which is significant only for very short chains^—. 
In the intermediate regime for 80 < N < 160 a Rouse-like 
behavior can be observed where the scaling exponent b is 
close to 1. For longer chains with N > 200 the value of 
b « 2.1 is typical of the crossover to reptation dynamics. 

The data in this plot are not closely following the re- 
ported scaling exponents because as the degree of poly- 
merization increases the density of the system changes 
(FigHJ). The different regimes are even less prominent in 
the MS MD simulations, where we observe a smoother 
transition between regimes than in UA MD data. In 
our calculations features like free volume effects due to 
chain ends (causing faster than Rouse decay of the dif- 
fusion coefficient for short chains) enter only indirectly 
as a consequence of the interplay of the input parame- 
ters, i.e. p, R gi N. The increasing of the density with N 
enters our analytical equation both directly, through the 
liquid density, and indirectly through the effective radius 
of gyration of the polymer, given that the compactness 
of the chain changes with the density (see also Section |Vl 
of this paper). 

The cis-l,4-PB chains are more flexible in comparison 
to PE chains, therefore we should expect larger entan- 
glement length N e (for PE from scattering N e ~ 130). 
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FIG. 4. Density as a function of degree of polymerization, N, 
for cis-l,4-Polybutadiene samples reported in Table |U 



Figure [3] shows that crossover from Rouse to reptation 
regime starts at N = 200. The largest system shown on 
Figure [3] is for N = 400 which is still a very weakly en- 
tangled system (1 or 2 entanglements per chain). Overall 
the agreement between simulated and reconstructed dif- 
fusion coefficients is quite good. 



IV. INTERNAL DYNAMICS 

The coarse-grained model adopted in this paper sim- 
plifies the macromolecular structure by reducing the de- 
scription of the molecule to an effective site that is cen- 
tered on the center-of-mass of the molecule. For this 
reason, the only dynamics that can be predicted by this 
CG description, through the rescaling procedure, is the 
long-time diffusive behavior of the center-of-mass of the 
molecule. However from the knowledge of the diffusion 
coefficient the monomer friction coefficient can be ob- 
tained, and the latter can be used as an input to the 
Langevin equation that describes the internal dynamics 
of the polymer chain at any length scale of interest. 

Specifically we follow our approach for the Cooperative 
Dynamics of a group of macromolecules, or Cooperative 
Dynamics Generalized Langevin Equation (CDGLE)^ 
Conventional theories of polymer dynamics, such as the 
Rouse theory for unentangled polymer dynamics and the 
"reptation" model for entangled dynamics, predict diffu- 
sive motion at short time following the ballistic regime; 
this is however in disagreement with simulations and ex- 
periments, which show a sub-diffusive regime before the 
com starts to follow Brownian motion. Both theories 
are mean-field approaches of single chain dynamics in an 
uniform bath. 18 Our approach focuses on the dynamics 
of a group of polymer chains in a melt and relates the 
anomalous subdiffusive behavior to the presence of coop- 



erative motion of a group of polymer chains in a dynam- 
ically heterogeneous liquid, as observed in simulations 
data of polymer melts and experiments^— Theoreti- 
cal predictions of the com and monomer mean-square 
displacements are shown to be in quantitative agree- 
ment with computer simulations of unent angle d 20 1 39 and 
slightly entangled polymer melts, 40 and with scattering 
experiments of Neutron Spin Echo. 21 

The physical picture underlying the theory builds on 
the fact that in polymer melts the dynamics appears het- 
erogeneous with the motion of a tagged chain correlated 
to the dynamics of a group of n chains comprised in- 
side the range of the intermolecular potential. The latter 
has a range of the order of the correlation hole, i.e. of 
the radius-of-gyration of the polymer, as it emerges from 
the solution of the OZ equation. The number of chains 
comprised in a volume of the order of is given by 
n « pTV 1 / 2 / 3 , where /, the statistical segment length, and 
the other quantities have been defined earlier on. The 
number of interpenetrating chains, n, increases with in- 
creasing density, degree of polymerization, and the stiff- 
ness of the polymer. 

In the Cooperative Dynamic approach the dynam- 
ics is described by a set of coupled equations of mo- 
tion (eom). Each equation is expressed in the space 
coordinates of the monomer a, belonging to molecule 
i and positioned at Ta\i), and contains a balance of 
different forces acting on the monomer. These in- 
clude the viscous force, C ~~j^~^ •> ^ ne intramolecular force 
— k s YlbLi Aa^r^ (t), which contains the structural ma- 
trix A, the time-dependent intermolecular mean- field 
force /r 1 i -^ln [rK<,ff(r (i) (<), r<*>(t))] , and the ran- 
dom interactions with the surrounding liquid, given by 
the random force Fa\t). 



n 



6=1 

1 d 



+F«(*) . 



In 



(14) 



Here A is the matrix of intramolecular connectivity, 
which reduces to the Rouse matrix when infinitely long 
and completely flexible macromolecules are considered, 
and is defined as 



N 



k,p=2 



(15) 



Here U is the equilibrium averaged bond correlation ma- 
trix 



< lfc • lp > 



(16) 
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and M is a structural matrix, with all the elements equal 
zero except = 1/TVfor i = 1,...,TV, = 1 with 
i = 2,..., TV, and Mi_ M = -1 with i = 2,...,' TV. The U 
matrix is a function of the local semiflexibility parameter 
g =< li • li+i > /(Wi+i), which is related to the persis- 
tence length of the polymer. For our samples the values 
of g are reported in Table [I] and are calculated from the 
molecule radius of gyration, which can be obtained from 
experiments or from atomistic simulations. 6 

Through Eq. ([T6]) the eom includes a complete micro- 
scopic description of the structure and local flexibility of 
the specific molecule under investigation, containing all 
the relevant parameters that define the intramolecular 
mean-force potential. Equations of motion for different 
beads belonging to the same chain or to two chains un- 
dergoing slow cooperative dynamics are coupled by the 
presence of intra- and inter-molecular interaction poten- 
tials. The intrinsic, chemical dependent semiflexibility of 
the macromolecule enters explicitly through the descrip- 
tion of the A matrix. 

The intermolecular potential is time dependent, as 
it is a function of the relative position coordinates of 
the center s-of-mass of a pair of molecules. As the two 
molecules move with respect to each other the inter- 
molecular force changes. From the initial ensemble of n 
dynamically correlated chains, molecules diffuse in time 
and finally move outside the range of the potential, R g , 
in a timescale of i 'decor -r = Rg/D. 

The set of coupled equations is solved after transfor- 
mation into normal modes of motion and numerically 
using a self-consistent procedure that calculate the ef- 
fective potential at any given distance. Once n and 
Cm = ksT / (ND) are defined, the first optimized against 
the simulation data and the second from our rescaling 
procedure, the set of equations has no adjustable param- 
eters. 

By applying consideration of symmetry the set of cou- 
pled eoms in normal coordinates can be simplified and 
reduced to two sets of TV uncoupled equations in the rela- 
tive and collective normal-mode coordinates . We assume 
that the eoms for internal modes (p > 0) do not contain 
intermolecular forces. This approximation is justified on 
the basis that polymer local dynamics is affected in a 
similar way by inter- and intra-molecular excluded vol- 
ume interactions, which in polymer liquids tend to com- 
pensate each other. Intermolecular forces, which enter 
the dynamics through the eom for the first (p = 0) nor- 
mal mode, still perturb the dynamics on the local scale 
through the linear combination of the normal modes. 

In the following, we present an overview of some of 
the quantities analyzed in the original simulations 15 and 
their comparison with the predictions of the Coopera- 
tive Dynamics theory having as an input parameter the 
reconstructed friction coefficient calculated from the soft- 
sphere simulations and the rescaling procedure as de- 
scribed in the first part of this paper. 

Figure [5] presents the end-to-end vector, time decorre- 
lation function for the TV = 112 sample of PB. The Coop- 



T 




t[ns] 

FIG. 5. End-to-end vector time decorrelation function for cis- 
1,4-Polybutadiene with N=112. Data from UA MD simula- 
tions (symbols) are compared with predictions of the theory 
for Cooperative Dynamics (solid line) where the number of 
correlated chains is set to n = 7 and the monomer friction 
coefficient is reconstructed from MS MD simulations, using 
the procedure described in this paper. 



erative Dynamic theory with the reconstructed monomer 
friction coefficient is directly compared against UA MD 
data. This function represents the rotational relaxation 
of the chain and cannot be measured from the soft col- 
loid simulation directly. The TV = 112 sample has chains 
shorter than the entanglement length and the Coopera- 
tive Dynamics theory reproduces the atomistic simula- 
tion data rather well. 

Figure [6] shows the center of mass mean square dis- 
placement (com MSD) for three samples of cis-1,4- 
Polybutadiene melts with chains of increasing degree of 
polymerization, TV = 240, 320 and 400. The long time 
diffusive dynamics is calculated from the diffusion coef- 
ficients obtained from the rescaling procedure. For time 
shorter than the longest relaxation time, r ~ R 2 g / D, the 
mean-square displacement of the com shows a subdiffu- 
sive behavior that cannot be reproduced by the Rouse 
approach, i.e. the conventional theory of unentangled 
chain dynamics. In the cooperative dynamics theory the 
subdiffusive dynamics is a consequence of the cooperative 
motion of the interpenetrating chains inside the correla- 
tion hole region, coupled by the effective intermolecular 
potential between the corns of the macromolecules. The 
theory, with the monomer friction from rescaled MS MD 
simulations as an input, is compared against the UA MD 
simulation data. The diffusion coefficient for PB240 is 
underestimated and for PB320 is slightly overestimated 
by the theory in comparison with the simulations, which 
is also can be seen in Table [TT] and Figure [3j These de- 
viations are within the error related to R g . Dot-dash 
lines in the inset of Figure [6] represent results calculated 
with the upper and the lower values of R g taken from 
the data of the atomistic simulations, as reported in Ta- 
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FIG. 6. Center of mass mean square displacement as a func- 
tion of time for cis-l,4-Polybutadiene samples with N=240 
(circle), 320 (square), and 400 (triangle). Predictions of the 
theory for cooperative dynamics (solid lines) are compared 
against the UA MD simulations (symbols). Also shown are 
the purely diffusive slopes obtained from the rescaled MS 
MD simulations (dashed lines). Inset illustrates how the un- 
certainty in the radius-of-gyration affects the mean-square- 
displacement for the N=240 sample: the upper and lower val- 
ues, reported as dot-dashed lines, corresponds to the upper 
and lower values of R g as reported in Table [I] 



ble U for the PB240 sample. The number of correlated 
chains is optimized to reproduce the subdiffusive behav- 
ior giving n = 10, 10, 12 for PB240, PB320 and PB400 
respectively, with a trend that qualitatively agree with 
the predicted scaling behavior. The error in the radius- 
of-gyration does not allow for a more precise calculation 
of the parameter n. 

Starting from the theory of cooperative dynamics, the 
single chain dynamic structure factor can be calculated 
as 

1 N 



i,3 
N 



1 E 

N 4^ 



exp 



.S_(( ri (*)-r,-(0)) 2 



.(17) 



EqfTTl describes the dynamics of a "tagged" polymer 
chain, which moves in a group of n interpenetrating 
chains. The intermolecular interaction between chains 
affects the dynamics resulting in the subdiffusive motion 
of the single-chain center-of-mass, as depicted in Figure 
[H and in the subdiffusive decay of the low q dynamics of 
S(q,t). 

The average time-dependent distance between two 
monomers in the tagged chain depends on the presence 
of the n — 1 other chains dynamically correlated, and it 
is expressed as a function of relative and collective mode 
coordinates. Once Qff and n are determined, the theory 
predicts the full decay of the dynamic structure factor, at 



FIG. 7. Normalized dynamic structure factor for cis-1,4- 
Polybutadiene with N=96 and q=0.04 (squares), 0.1 (circles), 
0.2 (triangles), 0.3 (diamonds) A -1 . The data from UA MD 
simulations (symbols) are compared against the Cooperative 
Dynamics theory (solid lines) where the number of correlated 
chains is set to n — 15 and the monomer friction coefficient is 
reconstructed from MS MD simulations, using the procedure 
described in this paper. 



the different experimental wave vectors. Figures [7] and 
[8] display the normalized dynamic structure factor ob- 
tained from the theory and compared with UA MD sim- 
ulations for q = 0.04,0.1,0.2,0.3 A -1 . For unentangled 
system PB96 shown on Figure [71 the agreement is excel- 
lent for all values of scattering vector q and time range. 
The Rouse model (not shown in this figure, see refJ^) 
fails to describe these data with satisfactory agreement, 
while the cooperative dynamics theory is found to be in 
quantitative agreement for samples at different molecu- 
lar weights, over the whole time range, and for different 
wave vectors. 

On Figure [8] the normalized dynamic structure factor 
is shown for the weakly entangled system PB400. The 
agreement for the smallest scattering vector which repre- 
sents global dynamics is excellent at all times, while for 
intermediate and large values of q the agreement between 
theory and simulations is slightly less accurate. The the- 
oretical decay, which is slightly faster than simulations 
at long times, suffers from the lack of properly describ- 
ing the entanglement effects. A more recent version of 
the cooperative dynamics theory includes chain uncross- 
ability, i.e. entanglements, and shows the correct slowing 
down of the relaxation in the long-time regime for entan- 
gled samples. 41 

Figure [9] displays the monomer mean square displace- 
ment for cis-l,4-Polybutadiene sample with N = 112. 
Theoretical prediction is compared against UA MD data 
from ref^. The agreement in both subdiffusive and dif- 
fusive regions is very good. The theory accounts for the 
local semiflexibility and chain connectivity, together with 
the cooperative motion of the chains due to the presence 
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FIG. 8. Normalized dynamic structure factor for cis-1,4- 
Polybutadiene with N=400 and q=0.04 (squares), 0.1 (cir- 
cles), 0.2 (triangles), 0.3 (diamonds) A -1 . The data from UA 
MD simulations (symbols) are compared against the theory 
for Cooperative Dynamics (solid lines) where the number of 
correlated chains is set to n — 12 and the monomer friction 
coefficient is reconstructed from the MS MD simulations, us- 
ing the procedure described in this paper. 



FIG. 9. Monomer mean square displacement, averaged over 
the innermost chain segments, as a function of time for cis- 
1,4-Polybutadiene with N=112. Predictions from the theory 
of Cooperative Dynamics (solid line) are compared against 
UA MD simulations (circles). The slope obtained from the 
rescaled MS MD simulations is shown as well (dot-dashed 
line) . 



of inter molecular interactions. 



V. THEORETICAL PREDICTIONS USING THE FREELY 
ROTATING CHAIN MODEL 

The purpose of the rescaling procedure is to predict 
diffusion coefficients from the mesocale simulations prop- 
erly rescaled when the atomistic simulations or the ex- 
periments are not available. One key input quantity to 
the procedure is the molecular radius of gyration. In the 
previous sections we adopted the values of R g determined 
from the atomistic simulations. When there are no data 
available for the radius of gyration, neither from atom- 
istic simulations nor from experiments, it is convenient 
to use some simple model, like the freely rotating chain 
model to provide R gi as it captures the semiflexibility of 
the polymer chain. In this section we present a critical 
discussion of adopting a freely rotating chain model to 
describe the PB chain, using a constant flexibility pa- 
rameter g. 

For the cis-l,4-polybutadiene chain with TV carbon 
atoms along its backbone there are uch=ch = N/4 dou- 
ble bonds with l^b = 1.34 A, ncH 2 -CH 2 — N/4 — 1 sin- 
gle bonds with = 1.54 A and ncH 2 -CH — N/2 single 
bonds with l' h = 1.5 A. The average bond length over all 
TV — 1 bonds is calculated as 
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FIG. 10. Density dependence of the semiflexibility parame- 
ter g calculated from Ea. ()19|) using values of the radius-of- 
gyration R g measured in UA MD and reported in Table Hill 



which in the N — >• oo limit gives the value of l% ve = 1.47A. 

The radius of gyration for a freely rotating chain 
(FRC) model with semiflexibility parameter g is given 



c>2 N 1 f]ave ,2 i 1 



2g 



(l- 5 )2 N 



(19) 

Figure [10] shows the density dependence of the semiflex- 
ibility parameter g in the UA MD simulations. 

Fitting simultaneously all samples to R^ A data gives 
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TABLE III. Semiflexibility parameter g calculated from FRC 
expression. 



System (R 2 ete ) U A [k 2 ] g" (gg^gj 



PB 32 


270±10 


0.6237 


45±5 


0.6237 


PB 48 


440±10 


0.6391 


70±5 


0.6243 


PB 56 


530±10 


0.6458 


85±7 


0.6339 


PB 64 


600±15 


0.6406 


95±10 


0.6246 


PB 80 


790±30 


0.6531 


125±10 


0.6374 


PB 96 


950±20 


0.6516 


152±16 


0.6394 


PB 112 


1150±20 


0.6609 


184±15 


0.6490 


PB 128 


1335±30 


0.6643 


215±18 


0.6544 


PB 140 


1430±40 


0.6576 


234±18 


0.6523 


PB 160 


1640±60 


0.6578 


275±25 


0.6595 


PB 200 


2100±80 


0.6634 


340±20 


0.6551 


PB 240 


2480±100 


0.6581 


410±30 


0.6557 


PB 320 


3100±200 


0.6382 


576±30 


0.6695 


PB 400 


4138±100 


0.6568 


678±30 


0.6519 



a R e te = R u et i (g ave = 0.651); b R ete = V6R% A (g ave = 0.645) 
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N 

FIG. 11. Radius of gyration squared over degree of polymer- 
ization Rg/N as a function of TV. UA MD data (circles), with 
statistical error, are compared against data calculated with 
the Freely Rotating Chain model using an averaged semiflex- 
ibility parameter, g = 0.6564 (solid line). 



an average value of the semiflexibility parameter for all 
densities, g = 0.6564. Figure fTTl displays the radius of gy- 
ration squared over degree of polymerization R 2 g /N as a 
function of TV from united atom simulations and from the 
freely-rotating-chain model. The curvature at small TV is 
due to chain-end effects. Samples with small TV present 
the largest deviation between the two sets of data, while 
the constant flexibility hypothesis works best for the large 
TV samples. The deviation is due to the fact that the 
united atom simulations are performed in the NPT en- 
semble where increasing TV corresponds to an increase of 
the liquid density, as shown in Figure |U The hypothesis 
of a constant semiflexibility parameter works well for the 
long chains for which coarse-graining methods are most 
useful. 

With this in mind, we report in Figure [12] the dif- 
fusion coefficient calculated from the mesoscale simula- 
tions, properly rescaled, where now the input radius-of- 
gyration is taken from the freely rotating chain model. 
As it is expected the best agreement between UA simula- 
tions and theoretical predicted diffusion is for the samples 
with larger TV. The disagreement at small TV is related to 
the fact that the density is changing with increasing TV in 
that region of the plot, and the radius of gyration should 
be properly corrected for this difference. In conclusion, if 
an exact value of the radius of gyration as a function of 
thermodynamic conditions and degree of polymerization 
is known, the described procedure should be able to pro- 
vide an accurate estimate of the dynamics starting from 
the mesoscale simulations of the coarse-grained system, 
which are computationally very efficient. The figure also 
reports a number of samples for which neither experi- 
ments nor simulations have been performed (TV = 180, 
280, 360, and 440). The most relevant conclusion of this 



o 




FIG. 12. Diffusion coefficients reconstructed from MS MD 
simulations using the radius-of-gyration calculated with the 
Freely Rotating Chain model (open circles) are compared 
against UA MD data (filled squares). Predictions for new 
systems with N=180, 280, 360, and 440 are shown as well 
(open diamonds). 



part of our study is the fact that it is possible to make 
predictions for data that are not reported in the UA MD 
simulations, obtaining information for new systems for 
which neither simulations nor experiments are available. 



VI. CONCLUSIONS 

In a couple of recent papers we have presented an 
original approach to rescale the dynamics measured in 
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mesoscale simulations of coarse-grained systems and re- 
cover the realistic dynamics as measured in atomistic 
simulations performed in identical thermodynamic condi- 
tions. The relevance of this method is in its ability of pre- 
dicting diffusion coefficient, with good precision, directly 
from the simulation of a coarse-grained system, without 
the need of performing an atomistic simulation. Sim- 
ply put, the method does not need a calibration curve. 
While the previous papers focused on the coarse-grained 
dynamics of polyethylene samples, this is the first pa- 
per in which the method is applied to a polymer with a 
different monomeric structure. 

In this paper we use an "extreme" level of coarse- 
graining as the whole chain is represented as a soft col- 
loidal particle. The choice of this representation is mo- 
tivated by two reasons. The first is that a possible er- 
ror in the procedure would be clearly made evident in 
this level of coarse-graining, and we believe that the pro- 
cedure would produce an even better agreement if the 
macromolecule is coarse-grained at a lower level, e.g. as 
a collection of soft particles, or blobs . 42 i 43 The second 
is that for this representation all the physical quantities 
that enter our approach are analytically solved and the 
rescaling procedure is formally derived. 

The paper presents a comparison between theoret- 
ical predictions of coarse-grained simulations properly 
rescaled following our first-principles procedure, and data 
from united- atom simulations for cis-l,4-Polybutadiene 
melts. The simulations were performed by Tsolou et al. 
in a NPT ensemble^ The rescaled diffusion coefficient 
shows good agreement with the diffusion measured in 
atomistic simulations, which supports the validity of the 
proposed procedure. 

As the coarse-graining model represents the whole 
polymeric chain as a soft colloidal particle, the only 
dynamical information that can be collected from the 
rescaled mesoscale simulation is the diffusion coefficient 
of the center-of-mass of the polymer. From this, how- 
ever, the friction coefficient of the monomer is derived. 
This paper shows how the monomer friction coefficient is 
used as an input to traditional Langevin equations for the 
dynamics of polymeric liquids. Specifically we use here 
the cooperative dynamics model for unent angled chains, 
given that the samples in the UA simulations are in the 
regime of unent angled and slightly entangled dynamics. 
From the theory of cooperative dynamics, which repre- 
sents the cooperative motion of a group of interacting 
polymers in a dynamically heterogeneous liquid, all time 
correlation functions of interest can be calculated. Specif- 
ically we present here the comparison with the quanti- 
ties reported in the original paper by Tsolou et al and 
show that the theory, with the rescaled friction coeffi- 
cient, is able to reproduce quantitatively the dynamics 
observed in the end-to-end time decorrelation function, 
and in the dynamic structure factor for unentangled and 
slightly entangled chains, and from the global to the local 
intramolecular dynamics. For entangled chains the local 
dynamics, represented by the high q regime in the dy- 



namic structure factor and in the monomer mean-square- 
displacement, shows a slight departure from the theory 
as the cooperative dynamics approach adopted here does 
not yet include the effect of entanglements. 41 
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